An Upper Mass Limit on a Red Supergiant Progenitor for the 
Type Il-Plateau Supernova SN 2006my^ 

Douglas C. Leonard^, Avishay Gal-Yam^, Derek B. Fox^, P. B. Cameron^, Erik M. 
Johansson^, Adam L. Kraus^, David Le Mignant^, Marcos A. van Dam^ 

ABSTRACT 

We analyze two pre-supernova (SN) and three post-SN high-resolution im- 
ages of the site of the Type Il-Plateau supernova SN 2006my in an effort to 
either detect the progenitor star or to constrain its properties. Following image 
registration, we find that an isolated stellar object is not detected at the loca- 
tion of SN 2006my in either of the two pre-SN images. In the first, an /-band 
image obtained with the Wide-Field and Planetary Camera 2 on board the Hub- 
ble Space Telescope, the offset between the SN 2006my location and a detected 
source ("Source 1") is too large: > 0.08", which corresponds to a confidence level 
of non-association of 96% from our most liberal estimates of the transforma- 
tion and measurement uncertainties. In the second, a similarly obtained V^-band 
image, a source is detected ( "Source 2" ) that has overlap with the SN 2006my lo- 
cation but is definitively an extended object. Through artificial star tests carried 
out on the precise location of SN 2006my in the images, we derive a 3 a upper 
bound on the luminosity of a red supergiant that could have remained undetected 
in our pre-SN images of logL/L© = 5.10, which translates to an upper bound on 
such a star's initial mass of 15 Mq from the STARS stellar evolutionary models. 
Although considered unlikely, we can not rule out the possibility that part of the 
light comprising Source 1 , which exhibits a slight extension relative to other point 
sources in the image, or part of the light contributing to the extended Source 2, 
may be due to the progenitor of SN 2006my. Only additional, high-resolution 
observations of the site taken after SN 2006my has faded beyond detection can 
confirm or reject these possibilities. 
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Subject headings: supernovae: general — supernovae: individual (SN 2006my) 
— stars: evolution 



Introduction 



The most common class of core-collapse supernovae (SNe) displays a distinct plateau in 



its optical light curve, and is therefore dubbed Type Il-Plateau (II-P; see lFilippenkdll997l for 
a review of SN classifications). This type of stellar explosion has long been thought to result 
from the core collapse and subsequent envelope ejection of isolated red supergiant (RSG) 



connection has begun to accumulate ( 


Van Dvk et al. 


2003; 


Smartt et al. 


2004; 


Maund et al. 


2005; 


Li et al. 


2006; 


Eldridee et al. 


2007; 


Schawinski et al. 


2008). 



By registering pre-SN and post-SN images, usually taken at high resolution using either 
space-based optical detectors, or ground-based infrared detectors equipped with laser guide 
star adaptive optics systems, progenitor star i dentifications have now been proposed for 
seven SNe II-P (for a contemporary review, see lSmartt et al.ll2008l . and references therein). 
Although different in detail, all seven proposed progenitors have properties consistent with 
those of supergiant stars. Because the field is still in its infancy — at this point none of the 
proposed progenitor objects has been definitively confirmed by its absence in high-resolution 
images of the SN site taken after the SN has faded beyond detection — it is imperative 
to carefully examine every new pro genitor c l aim. Here we investigate the status of the 
progenitor of SN 2006my, for which iLi et al.l (120071 . hereafter L07) recently proposed the 
identification of an RSG progenitor in pre-SN images and derived a zero-age main-sequence 
mass for it of Mzams = ^0- 



Although discovered several months after explosion (jStanishev &: NielsenI |2005| ). the 
photometry and spectroscopy of SN 2006my presented by L07 clearly establish it as an 
SN II-P. To identify the progenitor star, L07 register post-SN ground-based optical images 
taken with the Canada-France-Hawaii Telescope in the Sloan r' band under excellent seeing 
conditions (typical stellar full width at half-maximum [FWHM] = 0.6") with pre-SN Hubble 



^Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the Cahfornia Institute of Technology, the University of Cahfornia and 
the National Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. Additional data were obtained from the data archive at 
the Space Telescope Science Institute. STScI is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 
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Space Telescope {HST) Wide-Field and Planetary Camera 2 (WFPC2) images (typical stellar 
FWHM = 0.15" on the WF2 chip, on which the SN site is located), and identify a source in 
the pre-SN images within the 1 a error circle estimated from the transformation uncertainty. 
Here, we reexamine this identification with the benefit of two new sets of high-resolution, 
post-SN images: one taken with the wide-field channel of the Ne ar Infrared Camera 2 op er- 
ated behind the laser guide star assisted adaptive optics system ( IWizinowich et al.ll2006l ) on 
the Keck II 10-m telescope (stellar FWHM = 0.10"; note that SN 2006my served as its own 
tip-tilt star for the observations), and the other taken with the ifST WFPC2 camera (with 
the SN centered on the PC chip; typical stellar FWHM = 0.08") as part of a study of the 
progenitors of core-collapse supernovae (GO 10803; PI: Smartt). 

This paper is organized as follows. In § 12.11 we present the pre-SN and post-SN images 
and details of the photometric measurements performed on them; in § 12.21 we describe the 
image registration process along with our conclusion that no isolated stellar source is detected 
at the location of SN 2006my in either of the pre-SN images; in § 12.31 we estimate detection 
limits in the pre-SN images; and in § 12.41 we derive an upper mass limit on the progenitor of 
SN 2006my under the assumption that it was a RSG. We summarize our findings in § |3l 



2. Data Analysis 
2.1. Pre-SN Image Photometry with HSTphot 



Table [TJlists information on the five sets of image data considered by our study, hereafter 
referred to by the designation assigned to the final combined image from each data set. We 
first discuss our photometry of the two pre-SN images, VI and II. The individual frames 
comprising these datasets were preprocessed through the standard Space Telescope Science 
Institute pipeline using the latest calibrations as of 2008 July 23. The im ages were then 
further processed following the procedure described by [Leonard et al.l (120031 ). which employs 
the suite of programs desi gned specifically for the reduction of WFPC2 data that are available 
as part of the HSTphot (lDolphinll2000bl ) software package (ver. 1.1.7b; our implementation 
includes the most recent update of 2008 July 19). The hstphotta.sk automatically accounts 
for WFPC2 point-spread function (PSF) variations, charge-transfer effects across the chips, 
zeropoints, and aperture corrections. 

We performed photometry on the final, combined images using the hstphot task with 
option flag 10, which combines turning on local sky determination, turning off empirically 
determined aperture corrections (using default values instead), and turning on PSF residual 
determination. We ran hstphot on VI and II individually with a signal-to- noise ratio (S/N) 
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threshold of 1.0. 

In addition to flight-system magnitudes and uncertainties, hstphot returns several mea- 
surement parameters for each object detected. For our purposes, the most important of these 
are the "object type" and "sharpness" parameters . To determine the object type, hstphot 



compares the goodness-of-fit (i.e., the x value; see lDolphirul2000bl ) of the detected source's 
spatial profile with three different models: (1) the best stellar profile (i.e., the library PSF 
-|- residuals); (2) a single "hot" pixel without background; and (3) a completely fiat profile. 
If the stellar profile provides the best match, the object is labeled a "star" and designated 
as object Type '1', '2', or '3', depending on whether it is deemed a "good star", a "possible 
unresolved binary" , or a "bad star' (a star centered on a bad pixel), respectively. If the single 
"hot" pixel provides the best profile match, the object is labeled a "single-pixel cosmic ray 
or hot pixel" and designated as object Type '4'. Finally, if a flat profile matches best, the 



object is labeled an "extended object" and designated as object Type '5'. iDolphinI (j2000bl ) 
notes that because this is, by design, a conservative test for object type discrimination (i.e., 
it is a high threshold to be classified as anything other than a good star), many nonstellar 
objects will still be classified as Type 1 objects. For particular objects of interest, therefore, 
examination of the "sharpness" parameter is recommended, where sharpness values between 
-0.3 (object PSF broader than library PSF) an d +0.3 (object PSF narrower than library 



PSF) indicate confident point-source detections (IDolphinI l2000bl ). 



Figures [T]( a) and[T|^6) show an ~ 1" region surrounding the site of SN 2006my on pre-SN 
images II and VI (respectively), following transformation of both images to the pixel grid 
of image Kl (see § 12.21) . The SN 2006my site is rather complex, and while common objects 
are evident in Figures [T]( a) and[T]^6), the immediate vicinity of SN 2006my appears quite 
different in them. In II, an object of Type 1 ("good star") is reported by hstphot close to 
the SN location with a S/N of 5.6; we label this object "Source 1" in Figure [T^. In VI, an 
object of Type 5 ("extended object") is detected by hstphot with a S/N of 4.70; this extended 
source appears to occupy a "horseshoe-shaped" region that includes the SN location. We 
label this object "Source 2" in Figure [T]d. These two sources appear to be distinct from one 
another in the combined images. 

Source 1 is the object identified by L07 as the probable progenitor of SN 2006my. We 
confirm both the pixel location and photometry reported by L07 (also determined using 
HSTphot), with values measured by us (L07) of [x,y] = [410.23, 158.59] ([410.22, 158.63])! 



^Note that all pixel coordinates given in this paper are in the "IRAF" syste m — the coord inate system 
reported by using the imexamine task within IRAF or by running DAOPHOT IStetsonlligsrl ) — in which 
an integer value is assigned to a star that is centered on the center of a pixel. These pixel values are 0.5 
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and flight-system magnitude F814W = 24.48 ± 0.19 mag (24.47 ± 0.20 mag), which corre- 
sponds to / = 24.4 6 =b 0.19 mag (24.45 =fc 0.20 mag) fo l lowing transformation according to 
the prescriptions of iHoltzman et al.l (119951 ) and iDolphinI (j2000a|). As noted by L07, although 
formally classified as an object of Type 1 ("good star") by hstphot, Source 1 exhibits a slight 
east-west extension in the original WFPC2 image (see Fig. W^a); see also Figure 6 of L07). 
While L07 conclude that this object is most likely a single star, we note that hstpho t reports 
a "shar pness" value of —0.35 for it, which places it beyond the limits suggested by [Dolphin 
(l2000bl ) for confident point-source detections. Source 2 is unequivocally an extended source, 
as hstphot flags it as an object of Type 5 and measures its profile to have a sharpness of 
—0.58, well beyond the range for isolated sta r-like sources and for which hstphofs PSF star- 
fitting routines provide reliable photometry (jPolphin fc Kennicuttll2002l ). We note that this 
object may be the "extended source" mentioned by L07 but not further investigated due 
to its location near the ~ 2 a error radius of their derived SN pixel coordinates. We shall 
return to further discuss both of these "sources of interest" following a description of the 
image registration process. 



2.2. Image Registration 

To determine whether the progenitor star that exploded as SN 2006my is detected in the 
pre-SN images, we first registered image II to image Kl by using the IRAF tasks geomap and 



geotran, closely following the technique described by lGal-Yam et al.l (120051 ). To carry out the 
transformation, we used 16 common sources and obtained a final solution with a rms residual 
of 0.30 pixel in x and 0.27 pixel in y, in the Kl pixel grid. We then similarly registered images 
VI, V2, and 12 to the transformed II image (using more than a dozen common sources in 
each case), which resulted in five final images all registered to the common Kl pixel grid. 

Using the transformed images, we then measured the pixel locations of Source 1 (in 
image II) and SN 2006my (in images Kl, 12, and V2) using the centroiding algorithm of 
the imexamine task within IRAfI^I The results returned by imexamine depend somewhat on 
both the centering radius and fitting function employed. Thus, we applied a range of values 
— centering radii of from 3-10 pixels, inclusive, using both Gaussian and Moffat profiles 



greater in both x and y than those reported by hstphot, which foUows the DoPHOT (jSchechter et al.lll993f ) 



convention of assigning an integer value to a star that is centered in the lower left corner of a pixel. 

•^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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— and took the average as our "best" value and the measured scatter around the average 
as the measurement uncertainty. Due to the extended nature of Source 2 in image VI, it 
was not possible to estimate a precise pixel location for it by usi ng the IRAF centroiding 



algorithm; instead, we ran the source- finding program SExtractor (iBertin &: Arnoutslll996l ) 
on the image, which reported a nonstellar source and pixel position at the location of Source 
2. The pixel locations, measurement uncertainties, and transformation uncertainties of all 
objects are given in Table [2l 

Accounting for measurement and transformation uncertainties, we determine offsets 
between Source 1 in image II and SN 2006my in images Kl, 12, and V2, of 0.083" ± 0.017", 
0.080" ±0.016", and 0.086" ±0.016", which correspond to separations significant at the 4.8 o", 
4.9 0", and 5.2 a levels, respectively, for a two-dimensional (i.e., both x and y) Gaussian. 
This indicates non-association at greater than the 99% confidence level (see the Appendix 
for a full discussion of how significance of source separation is determined), implying that 
the object detected as Source 1 is in all likelihood not the progenitor of SN 2006my. Due to 
the extended nature of Source 2, it is not possible to derive a similarly well- quantified offset 
and uncertainty estimate between it and SN 2006my, but we note that from the location 
reported by SExtractor, its formal separation is only 0.024", and examination of Figure [T] 
reveals that the location of SN 2006my is indeed coincident with part of the extended region 
identified as Source 2. 

To serve as a check on the results obtained by transforming all images to the Kl frame, 
we also performed a direct, "HST-only", registration between images V2 and VI, and 12 and 
II, using the /t-stp/tot- reported pixel positions (with an additional correction for distortion 



using the solutions of lAnderson fc King] 120031 ) for all (> 20) objects used to determine the 



transformation. Results of these registrations are given in Table [31 To assign a meas urement 



uncerta inty on the /istp/iot-measured object pixel positions, we consulted Figure 4 of [Dolphin 



( l2000bl ). which presents estimates of the 1 a astrometry error for sources detected by hstphot 
as a function of count level. For the count level of Source 1 (~ 45), a total astrometry 
error estimate of ~ 0.4 WF2 pixel, or ~ 0.29 pixels in both x and y, is determined, while 
for SN 2006my (counts > 20, 000) position uncertainties of only 0.03 pixels in x and y are 
derived. 

Accounting for measurement and transformation uncertainties, we determine the sepa- 
ration between Source 1 in image II in the transformed location of SN 2006my from image 
12 to be 0.082" ± 0.031", which represents an offset of 2.6 a for a two-dimensional Gaussian, 
indicating non-association at the 96% confidence level (see the Appendix). This separation 
is nearly identical to those determined in the Kl pixel frame, although the significance of the 
offset is lower. This is due mainly to the larger uncertainty adopted for the pixel location 
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of Source 1 in II than was found empirically for the II image transformed to the Kl pixel 
grid. Registration of II onto Kl resamples the data onto a finer grid and therefore likely 
provides a better localization of Source 1, although, as noted by L07, resampling the data 
has the risk of smoothing several (extended) sources into a point source. We also note that 
a direct registration of V2 to II (not listed in Table 3) yields an offset between Source 1 
and SN 2006my of 0.092" ± 0.032", which represents an offset of 2.9 a for a two-dimensional 
Gaussian and indicates non-association at the 98% confidence level. With all of these lines 
of evidence, we conclude that a 96% confidence level for non-association is a conservative 
estimate of the significance of the offset between Source 1 and SN 2006my. 

We note that in the original II WFPC2 image. Source 1 presents itself as elevated 
flux in two pixels (i.e., pixels [411, 158] and [410, 159]; see Figure 6 of L07), one of which 
contains our derived location of SN 2006my at its far edge (Table [3]). The slight extension 
of Source 1 noted earlier (§ 12.11) raises the possibility that this source may, in fact, contain 
light from more than one object — for instance, two RSG lying in adjacent WF pixels, 
or perhaps a compact star cluster — an idea considered by L07 but deemed improbable. 
This lingering possibility can only be definitively removed from consideration by future 
high-resolution imaging after the SN has faded beyond detection. For now, we conclude 
that from our astrometric measurements, an isolated stellar object is not detected at the 
location of SN 2006my in either of the two pre-SN images: In image II, the offset between 
the SN 2006my location and the source detected by hstphot is too large, and in VI the 
source that overlaps with the SN 2006my site is extended. The relevant question that arises, 
then, is this: How faint a point source at the precise location of SN 2006my would we have 
confidently detected as a point source in the pre-SN images? 



2.3. Detection Limits in Pre-SN Images 

To set accurate detection limits on a point source at the location of SN 2006my in the 
pre-SN images, we proceeded as follows. First, we used the showpsf task within HSTphot to 
produce the library PSFs appropriate for point sources at the precise pixel locations (accurate 
to 0.1 pixel) of SN 2006my in the VI and II images. Using these PSFs, we then inserted 
artificial stars of known flux (corresponding to 21.5 < < 27.0 and 20.5 < I < 26.0) at the 
SN locations in the VI and II images. Then, we ran hstphot on these images in exactly the 
same manner as we did when seeking a progenitor star in the original images. 

To be considered a confident "detection" of a single star at the SN 2006my location, 
we demanded that the object be classified by hstphot as a "good star" (Type 1), have a 
sharpness between -0.3 and -1-0.3, and have a reported pixel location no more than 1 a 



- 8 - 



away from the nominal SN location, where the uncertainty accounts for the measurement 
and transformation uncertainties listed in Table [3] as well as the astrometric uncertainty 



appropriate for the object from Figure 4 of iDolphid (l2000bl ) 



From this analysis, we derive detection limits of = 24.9 ± 0.3 mag and / = 24.4 ± 0.2 
mag for point sources at the location of SN 2006my in the VI and II images, respectively. 
In both cases the limiting magnitude is set by the sharpness parameter becoming less than 
-0.3 (i.e., source too extended to be confidently considered a point source). We note that 
these detection limits are significantly shallower than those derived by L07, in which detection 
limits were derived by examining the magnitudes of all 3 o detections in the images. Because 
we are specifically interested in the point-source detection limits at the location of the SN, 
we consider our (less restrictive) limits to more accurately reflect the relevant detection 
threshold. 



2.4. Properties of the Progenitor Star 



To convert our detection thresholds to constraints on the initial mass of an RSG progeni- 
tor star that could have exploded as SN 2006my, we employ the metallicity-dependent stellar 
models produced with the Cambridge stella r evolution code, STARS, the descendant of the 
code d eveloped originally by EggletonI (Il97ll ) and updated most recently by lEldridge &: Tout 



( I2OO4J . see also ISmartt et al. 



20081 . and references therein) The models follow stellar evo 



lution up to the initiation of core neon burning, which is likely to give an accurate indi 
cation of the pre-SN l uminosity; comparisons with other contemporary model grids (i.e. 
Heger fc Langen I2OOOI : iMeynet fc Maederl |2000| ) show that the endpoints fo r stars in the 
8 ^ 15 Mq range differ by at most 0.2 dex in luminosity among the codes (jSmartt et al. 



20041 ). 



Using the metallicity and radial gradient in NGC 4651 published by iPilyugin et al. 



f l2004h . L07 derive a metalhcity at the SN 2006my loca tion of logfO/jj") + 1 2 = 8.51 ± 
0.06, which is subsolar according to the recent analysis of lAsplund et al.l (120041 ). who found 
[OjE] = 8.66 ± 0.05. We thus use the Z = 0.01 (the closest metalhcity calculated) STARS 
stellar evolution models as the basis for deriving the upper mass limit for our study. Figure [2] 
displays the final predicted luminosity for stars with masses between 8 and 20 Mq from these 
models. 



We now seek to determine the lowest possible luminosity that an RSG could have had 



"'The models were downloaded from the code's Web site, at |http: / /www.ast.cam.ac.uk/^stars 
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and still have been confidently detected by our analysis of the pre-SN images. Since it 
will prove to be far more restrictive, we begin by considering the II image, and calculate 
the bolometric magnitude of RSG stars corresponding to our detection limit through the 
equation: 

Mbol = -fi - Av + I + (V - J)rsg + BCy, 

where /j, is the distance modulus of NGC 4651, Ay the extinction to SN 2006my, I the J-band 
detection threshold, (V — /)rsg the color range of RSG stars (i.e., spectral types K3 M4), 
and BCv the bolometric correction corresponding to each (V — /)rsg- As detailed by L07, 
distance estimates to NGC 4651 span quite a wide range. Because it is our goal here to set 
the most conservative lower bound on ou r detection threshold , we adopt the long distance 



estimate, /i = 31.74 ± 0.25 mag, found by ISolanes et al.l (120021 ) by averaging seven different 



TuUy-Fisher distances to the galaxy. For the extinction, we adopt Ay = 0.08 ± 0.02 mag, 
which represents the Galactic value along the line of sight, because there is no evidence 
for host-galaxy extinction (L07). For the co lor and bolometric corrections appropriate for 



RSG stars, we consult the values reported by lElias et al.l (1 19851 ). who find that the quantity 



[{V — /)rsg + BCy] lies in the remarkably tight range 0.88 — 1.0 for supergiant stars of 
spectral types K3 - M4 (i.e., RSG). Again, because we wish to set the most conservative 
lower detection limit, we assign (V — /)rsg + BCy = 0.88, which is the value obtained for 
an M4 supergiant star. Finally, we set / = 24.4 ± 0.2 mag, as derived in § 12.31 

With these values, we derive Mboi = —6.54 ±0.32 mag as the limiting bolometric magni- 
tude, above which any RSG would have been detected in our pre-SN image. This corresponds 
to a 3 cr detection threshold of Mbd = —7.50 mag, which translates to logL/L© = 4.90. If 
we adopt a maximum systematic uncertainty of 0.2 dex in the theoretical stellar model end- 
points, then the final 3 a lower bound on the luminosity of a RSG that would have been 
confidently detected in our pre-SN image is logL/L0 = 5.10. From Figure [2] this corre- 
sponds to an upper bound on the progenitor mass of Mzams = 15 Mq, and we therefore 
conclude that any RSG progenitor with an initial mass greater than this value would have 
been detected using our analysis procedure. 

Applying a similar analysis to the VI image results in a detection threshold of log L / Lq = 
6.30, which unfortunately does not rule out any progenitors up to 200 Mq, the highest pro- 
genitor mass considered by the STARS models. Our most restrictive limit thus comes from 
the II image. 



3. Conclusions 



We analyze two pre-SN and three new post-SN high-resolution images of the site of the 
Type Il-Plateau supernova SN 2006my in an effort to either detect the progenitor star or 
to constrain its properties. Our primary result is that we do not detect an isolated stellar 
object at the location of SN 2006my in either of the two pre-SN images. Prom our image 
registration, we therefore do not confirm the association found by L07 between a stellar 
source (Source 1) in pre-SN 7-band images (II) and the location of SN 2006my. Using new 
high-resolution post-SN images, we derive an offset between SN 2006my and SoTirce 1 of 
> 0.08" from the SN location, which represents a confidence level of non- association of more 
than 96% from our most liberal estimates of the image transformation and measurement 
uncertainties. Through artificial star tests carried out at the precise location of SN 2006my 
in image II, we derive a 3 a upper bound on the mass of the progenitor of SN 2006my of 
-^ZAMS = 15 Mq from the STARS stellar evolutionary models. 

Although considered unlikely, we can not rule out the possibility that part of the light 
comprising Source 1, which exhibits some extension relative to other point sources in image 
II, or part of the light contributing to Source 2, a definitively extended source detected in 
pre-SN y-band images (VI) that has overlap with the SN 2006my location, may be due to 
the progenitor of SN 2006my. Only additional, high-resolution observations of the site taken 
after SN 2006my has faded beyond detection can confirm or reject these possibihties. 
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A. APPENDIX 

ON THE CONFIDENCE LEVEL OF NON-ASSOCIATION OF POINT 

SOURCES 



Determining the significance of positional offsets between objects identified in two images 
is of paramount importance when assessing the potential association between a progenitor 
star (in a pre-explosion image) and a supernova (in a post-explosion image), where one 
image has been transformed into the pixel frame of the other. Six sources of uncertainty are 
typically identified: Positional uncertainty of the putative progenitor star (in both x and y) , 
positional uncertainty of the SN (in both x and y), and the uncertainty in the transformation 
(also in both x and y). In much of the past work involving SN progenitors (hereafter, the 
"traditional" approach), offset significance (i.e., how many "sigma" away the two objects 
are) has been derived by taking the measured radial offset between the two objects, and 
then dividing this value by the quadrature sum of all of the uncertainties. A formal analysis, 
however, finds this approach to be somewhat in error. 

Consider the two-dimensional Gaussian that describes the uncertainty in an SN's mea- 
sured position {x,y) compared with that of a putative progenitor star {xo,yQ): 

p{x,y)dxdy = ^-{x-x,r/2ai^-(y-yo)V2al ^^^y^ 

where o^^y^ represents the quadrature sum of the measurement and transformation errors in 
x{%j). Converting this to a normalized form by defining E^, = (x — xo)/axi T^y = {y — y{))/(7yi 
yields: 



p(S^,E^)(iE^c?Sj, = _l_e-(s'/2)e-(4/2)(a^(iE^)((jj^o?S^), 
which becomes the expected: 



2nax(Ty 



p(E„ Ej,) ciE, dEy = i-e-(^'/2)e-(^?/2)dE, dEy. 

We convert this normalized Gaussian into radial form with a further change of coordinates 
(Ea, = pcos9, T,y = psinO) and a little algebra to yield: 



p{p, 6) dp de = ^e-^'/^p dp de, 
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or ignoring the angular part of the distribution, 

p{p) dp = pe^^^^^ dp. 

The total integrated probability for p < Q (where Q is equal to, e.g., [S^ + S^]^/^), then, is 

p(p<g) = l-e-(«'/2)_ ^A^) 

A simple example serves to highlight the differences between calculating offset signifi- 
cance in the "traditional" manner versus calculating it according to equation (1A1|) . Suppose 
a progenitor star candidate and a SN are measured to be located at pixel locations [10,10] and 
[11,11] on pre-SN and (transformed) post-SN images, respectively, with total uncertainties 
(quadrature sum of all measurement and transformation uncertainties) of = 1, Cy = 1. In 
the traditional approach, this would represent an offset of 1.41 ± 1.41 pixels, or a separation 
significant at the la level, which implies a confidence level of non-association of 68%. 

However, if we calculate the total integrated probability for p < 1.41 according to 
equation flAip we find p{p < 1.41) = 1 — exp[— (1.41^/2)] = 0.63, which implies a confidence 
level of non- association of 63%. Because we are dealing with positions and uncertainties in 
two dimensions, it thus appears most precise to state the result as follows: "The two objects 
are located lAla away from each other for a two-dimensional Gaussian, which indicates 
non-association at the 63% confidence level." We therefore express all of our association 
significances in this paper in this manner. 
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Table 1. Observations of NGC 4651 



Telescope 


Instrument 


Data Archive Name 


UT Date 


Exposure Times 
(s) 


Filter ( 


Plate Scale 
arcsce pixel ^) 


FWHM 
(arcsec) 


Combined Image 
Designation 


HST 


WFPC2 


u2dt0901/2/3t 


1994 May 20 


60 + 300 + 300 


F555W 


0.10 


0.15 


VI 


HST 


WFPC2 


u2dt0904/5/6t 


1994 May 20 


60 + 300 + 300 


F814W 


0.10 


0.15 


11 


Keck II 


LGS/NIRC2 




2006 Nov 28 


10 X 6 X 10 


Kp 


0.04 


0.10 


Kl 


HST 


WFPC2 


u9ox0301/2/3/4m 


2007 Apr 26/27 


4 X 300 


F555W 


0.05 


0.08 


V2 


HST 


WFPC2 


u0ox0305/6m 


2007 Apr 27 


500 + 700 


F814W 


0.05 


0.08 


12 


Note. — 


SN 2006my in NGC 4651 was discovered 


on 2006 November i 


3.82 UT jNakano 


& ItaKakill200eD 


. Pre-explosion 


HST/WFPC2 


images {VI and 11) were 



obtained as part of a kinematic study of the core of NGC 4651 {GO 5375; PI: Rubin). Post-explosion image Kl was reduced according to the methods detailed by 
I Gal- Yam et al ] i2005lV Details on the reduction and analysis of all images obtained with HST/WFPC2 are given in the text. Plate scale is pixel size of the detector 
chip on which the site of SN 2006my is located {i.e., WF2 for VI and II, PC for V2 and 12, and wide-field channel for Kl), and FWHM is the measured full width 
at half-maximum of a point source in each final, combined image. We note that the scale of the wide camera (40" X 40" ) is larger than the isoplanatic patch size, 
which yields a variable PSF shape over the full image. 
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Table 2. Image Transformations to LGS (Image Kl) Pixel Grid 



Measurement Uncertainty Transformation Uncertainty 



Image 


Object 


X 


Y 








ay 


11 


Source 1 


579.55 


543.28 


0.34 


0.04 






Kl 


SN 2006my 


581.55 


542.72 


0.05 


0.03 


0.30 


0.27 


12 


SN 2006my 


581.47 


543.07 


0.05 


0.02 


0.21 


0.22 


V2 


SN 2006my 


581.69 


543.07 


0.02 


0.02 


0.23 


0.25 


VI 


Source 2 


581.95 


542.55 






0.15 


0.15 


Note. 


— Location of 


sources 


identified 


in Fig 


. [1] in images 


listed in Tabled! following 


; transformation of 



all images to the Kl pixel grid. All transformation uncertainties are given in pixel units (with a plate scale 
of 0.04" pixel"^; see Tabled]) and represent the rms residuals following registration of each image to image 
II transformed onto the Kl grid. Because Source 2 is extended (and not symmetric), we do not assign a 
measurement uncertainty on its position. Measurement uncertainties were derived as discussed in the text 
(§ 12. 2p . Note the particularly large uncertainty in the x coordinate of Source 1 in transformed image II, 
likely due to the slight extension of the object along the x (i.e., east- west) direction, as discussed in the text 
(§ [12]) and evident in Fig. dJi- 



Table 3. Image Transformations to 1994 WFPC2 (Images VI and II) Pixel Grid 



Measurement Uncertainty Transformation Uncertainty 



Image 


Object 


X 


Y 


fx 


(Jy 


fx 


ay 


11 


Source 1 


410.23 


158.59 


0.29 


0.29 






VI 


Source 2 


410.31 


159.58 










12 


SN 2006my 


410.18 


159.41 


0.03 


0.03 


0.15 


0.12 


V2 


SN 2006my 


410.20 


159.51 


0.03 


0.03 


0.12 


0.09 



Note. — Location of sources identified in Fig. [T] in images defined in Table 1. Transformation uncertain- 
ties are given in pixel units (with a plate scale of 0.10" pixel"^; see Table [1]) and represent the rms residuals 
following registration of 12 onto II and V2 onto VI. The locations of Source 1 and Source 2 in images II 
and VI, respectively, are those reported by hstphot. The locations of SN 2006my in the VI and II pixel 
grids were derived by taking the hstphot-ieported coordinates of SN 2006my in images V2 and 12 and trans- 
forming them using the geoxytran task in IRAF to the VI and II pixel gr ids, respectively. Measurement 
uncertainties for all object locations are those suggested bv iDolphinI (j2000bl ) . Because Source 2 is extended 
(and not symmetric), we do not assign a measurement uncertainty on its position. 
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Fig. 1. — Site of SN 2006my in pre-SN images II (panel a; all image designations are as 
given in Tabled]) and VI (panel b), and post-SN image 12 (panel c). All images have been 
transformed and resampled to the pixel grid of the Kl image (§ l2.2p . The cyan circles indicate 
the approximate 5a uncertainty (for a two-dimensional Gaussian; see the Appendix) of the 
position of SN 2006my relative to the transformed II image as measured in the transformed 12 
image. The white circles indicate the same level of uncertainty in the location of SN 2006my 
as measured in the Kl image relative to the transformed II image. Two "sources of interest," 
discussed in the text, are labeled Source 1 and Source 2 in panels a and b, respectively. 
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Fig. 2. — Initial mass vs. final predicted lumino sity prior to explosion for Z = 0.01 stars 
evolved with the STARS stellar evolution code (jEldridge &: ToutI l2004l ). The dashed line 
indicates the 3a upper luminosity limit for a RSG that could have remained undetected by 
our analysis of pre-SN images of the site of SN 2006my. 



